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We present a microscopic theory of bipolar quantum well structures in the photovoltaic regime, 
based on the non-equilibrium Green's function formalism for a multi band tight binding Hamiltonian. 
The quantum kinetic equations for the single particle Green's functions of electrons and holes are self- 
consistently coupled to Poisson's equation, including inter-carrier scattering on the Hartree level. 
Relaxation and broadening mechanisms are considered by the inclusion of acoustic and optical 
electron-phonon interaction in a self consistent Born approximation of the scattering self energies. 
Photogeneration of carriers is described on the same level in terms of a self energy derived from 
the standard dipole approximation of the electron-photon interaction. Results from a simple two 
band model are shown for the local density of states, spectral response, current spectrum, and 
current-voltage characteristics for generic single quantum well systems. 
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I. INTRODUCTION 

Since the pioneering work of Barnham and co-workers^ 
in the early nineties, the potential efficiency enhancement 
by the introduction of quantum wells in the intrinsic re- 
gion of a pm-diode solar cell (Fig. [1]) has attracted con- 
siderable interest both from the photovoltaic community 
and within a broad spectrum of fundamental research 2 . 
A consistent and quantitative description of the carrier 
generation, recombination, relaxation and transport pro- 
cesses in quantum well solar cells (QWSC) requires the 
combination of a microscopic model for the electronic 
structure with a formalism for quantum transport in in- 
teracting systems. The non-equilibrium Green's func- 
tion formalism (NEGF), first introduced by KadanofF and 
Bayrn^ and by KeldysbJ, together with a tight-binding 
or Wannier basis meets these requirements and has been 
succssfully applied to similar systems such as quantum 
cascade lasers^, infrared photodetectors 7 -, carbon nan- 
otube photodiodes&2 or resonant tunneling diodes^. 

The paper is organized as follows: In section II, we 
introduce the model Hamiltonian in a planar orbital ba- 
sis and the procedure based on the NEGF formalism to 
use it in the derivation of physical quantities. Section 
III presents and discusses typical results of the theory 
for a generic bipolar quantum well structure. Section V 
summarizes the paper and provides an outlook to future 
work. 



II. MICROSCOPIC MODEL FOR QWSC 

A. Hamiltonian and basis 

The QWSC system is described in terms of the model 
Hamiltonian 



H — Hq + Hi, Hi — H, 



ep ~r H e j . 



(1) 




Hq provides ballistic transport: it contains the kinetic 
energy, the (bulk) band structure and band offsets, and 
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FIG. 1: (color online) Characterizing structure and processes 
of a pin-QWSC. 



also includes the electrostatic potential from the solution 
of Poisson's equation, which corresponds to the consid- 
eration of single species carrier-carrier scattering on the 
Hartree level. The interaction part Hi consists of the 
terms H ep and H ei for electron-phonon and electron- 
photon scattering, respectively. Any other kind of inter- 
action, like scattering by ionized impurities, alloy com- 
position inhomogencitics or interface roughness, inter- 
carricr-scattering beyond the Hartree level, and all non- 
radiative recombination processes like Auger or trap re- 
combination are contained in additional terms that will 
not be discussed in this paper, which is focused on the 
radiative limit. 



In layered semiconductors, carrier Bloch states can be 
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represented in terms of linear combinations of planar 
orbitalsii 

|n,(k,fc,)> =5^C a , £ (k J fe,)|a,L,k), (2) 
|a,L,k) = ^=5>* k ^|a,L,Rf), (3) 

R H 

where n is the band index, a denotes a set of orthogonal 
localized orbitals (e.g. s,p x ,p y ,p z , s* in a 10 band model 
for zinc-blende materials^), L indicates the layer, which 
can consist of several different atomic layers, and Rj^ the 
(transverse) location within the layer. N is a normaliza- 
tion factor and k, k z the transverse and longitudinal wave 
vectors, respectively. The corresponding field operators 



k,L a 

A klr)^ 

k,L a 



,L.k> 



(4) 
(5) 



where Ca,i,k (c„ ^ k ) is the annihilation (creation) oper- 
ator for a fermion in state |a,L,k). 

In a planar orbital basis (POB), the Hamiltonian for 
ballistic transport is expressed as 

^ = EEE (k)(l - 5i',i)4,L,k 6 «',i',k 



k a, a' L,L' 
Da,L;a',L(k)5 



L,L>C aLV C a >,LM\ 



(6) 



where D contains the on-site energy, the intra-layer cou- 
plings (overlap integrals) and the Hartree potential, while 
t denotes the inter-layer coupling. 

The operator for carrier-photon interaction reads 



V = —A ■ p, 

m 



(7) 



with the quantized photon field given by 

&A,q(*) 



1 



Aq 



eA q e 



vqr 



(8) 

where eA q is the polarization of the photon in mode A 
and with momentum q created by the boson creation 
and annihilation operators {&',&}, and V is the absorbing 
volume. 

In a first approach, we restrict the discussion to single- 
mode monochromatic photons of energy hujj and use the 
standard dipole approximation, which yields^ 



A = 



I Hy/fIe<f) Uy 



N 7 c 



(9) 



where a is the polarization and ^ represents the in- 
coming photon flux, which depends on the intensity J 7 
and the photon energy, and provides iV 7 photons per ab- 
sorbing volume V and for given optical properties (e: di- 
electric constant, /x: magnetic permeability) . In the POB 
for a layered system, the Hamiltonian for electron-photon 
interaction takes the form 



H ei = J d?r^(r)V ei ^{v) (11) 
= EEEKw( k )lL,k^.i'* 



L.L' a,ot' k 



x (be 



-lUJ-yt 



(12) 



In the dipole approximation, A has no spatial depen- 
dence and thus 

MZl-cc'L'W = — A (a,L,k|p|a',L',k), (13) 



TO 



-a, e is the electron charge and 



where An = , 

too its bare mass. The band structure model depen- 
dent dipole-matrix elements for the (direct) interband 
transitions can be written in terms of the tight-binding 
Hamiltonian a o 13 ' 14 i 15 

1 



{a,L,k\p\a',L',k) 



E 

Rf,Rf 



»k||-(R|, -Rm ) 



x (a,L,Rj^|p|a',L',Rf 



(14) 



(a,L,R[|p|«',L',R[') =^%,X,R 



r,Hr 



= — 7-(R L — R L )[-ffo]a.L;a'.L' 

in 



(15) 

where R L = (R^, LA). In the case of light incidence 
normal to the layer, the polarization is purely transverse, 
and M 7 becomes a scalar function of the transverse mo- 
mentum. 

For the interaction of carriers with phonons, which is 
on the level of a coupling to an equilibrium heat bath, the 
harmonic approximation provides the interaction term 



(16) 



where U q characterizes the coupling matrix elements. In 
the case of a diatomic basis, as in zinc-blende compounds, 
the corresponding POB interaction Hamiltonian is given 
by 



H ep = / d 3 rft(r)V ep $(r) 



(17) 

= E E E M a P L(<l)cl,L,kCa,L,k- 9( 
L,k q a 

x (6q + a- q )- (18) 

where the exact form of the coupling element M ep de- 
pends again on the band structure model. 
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B. Green's functions and self-energies 

Within the planar orbital basis, the real time non- 
equilibrium Green's functions are defined as the non- 
equilibrium ensemble averages 

G< L;a ,, L ,(k;t,t') = J(4,i',k(* , )^,i,k(*)> J (19) 

GlwAW) = -l(^.L.^t)cl Llk (t')), (20) 
Glw (k; t, t>) = 0(t - H)[G> L;a , i£ , (k; t, t') 

-G< L , al ^{W)l (21) 
Gi, Lv *, v (k; t,t') = 9(t' - i) [G< L;Q ,, L , (k; t, t') 

-G> iL . a , tL ,(k;t,t')}. (22) 

In steady state, the above Green's functions depend only 
on the time difference r = t — t' , and it is thus possible 
to work with the Fourier transform 

G a>L . a ,, v {k-E) = I ' dre iET ' h G a , L , al , L ,{\<i;r), r = t-t'. 

(23) 

The effects of carrier injection and absorption by ex- 
tended, highly doped contacts acting as reservoirs are ab- 
sorbed into respective boundary self energies S B , reflect- 
ing the openness of the system and leading to an effec- 
tive Hamiltonian of the truncated syste m 16 i 17 i 18 . Since 
the contacts form equilibrated flat band regions, their 
propagating and evanescent bulk Bloch states can be de- 
termined exactly. The boundary self energy then rep- 
resents the matching of the planar orbital states in the 
device to the extended lead modes at the interface of 
the contacts, corresponding to a quantum transmitting 
boundary method^. For instance at the left boundary 
(L = 1), the retarded boundary self energy is given by 
(see appendix for a detailed derivation) 



Z?f(k,E) =t 1;0 (k) 



(U^,E)[A z (k,E)]- l [U.(^,E)}- 



(24) 



where U— specifies the transformation from localized ba- 
sis to left-travelling Bloch states and A z is the interlayer 
propagator for the corresponding bulk mode a 20 ' 21 i 22 . The 
lesser and greater self energies are then obtained from the 
broadening function Tf and the Fermi distribution j jlL 
of the contact characterized by the chemical potential 



£<f(k,£) =tf, L (E)r?(k,E), 
Z>?(k,E)=-i[l-U L (E)]T?(k,E), 

rf(k,£) =i\E**-(i? 



(25) 
(26) 
(27) 



Analogous expressions arc found for the right contact. 

While the boundary self energies result from an ex- 
act treatment, interactions such as carrier-phonon and 
carrier-photon scattering are included pcrturbatively in 



terms of interaction self energies £ on the level of a self- 
consistent Born approximation (SCBA). The self ener- 
gies for both carrier-photon and carrier-phonon are ob- 
tained from the Fock term in second order perturbation 
theory for general carrier-boson interaction. The corre- 
sponding Hartree term is neglected at the present stage 
(see e.g. [23| for an extensive discussion). In the case of 
the light-matter interaction the Bosc-Einstcin distribu- 
tion Nq^hujq) in the equilibrium bosonic propagator is 
replaced by the number of photons iV 7 present in a layer. 
The lesser and greater self energies read (in full matrix 
notation) 

Sf 7 (k; E) = ihAP(k) [/V 7 G^(k; E^hw^ 

+ (iV 7 + l)G^(k;£±/iw 7 )]M 7 (k), (28) 
and the retarded self energy is given by 



£f 7 (k;£) =jfiM^(k) 



+ N 1 G R {k;E + hw 1 ) + -[G<(k;E- Hlu 7 ) 
-G < (k;E + huJ 1 )} 

dE' ( G<(k;E-E') 




2tt 



G<(k;E-E') 
E' + fiw 7 



E> - tvjj~ 



AP{k) 



(29) 



The principal value V in the expression for the retarded 
self energy is often neglected, since it will only con- 
tribute an energy renormalization, but not to relaxation 
or phase breaking. We will adopt this approximation in 
the present work. 

For the interactions with polar optical phonons, the self 
energies are then given by (again neglecting the principal 
value integration in the retarded case) 

^S, L ' (k; £)=E MPOp ( k > Hi i ^ «! L '> «') 



[N LO Gf L . a , LI {ci\\-ET^Lo) 

(N LO + l)Gt !L . a , !L ,(c ir ,E±f^ L o)] 

(30) 

M pop {k,^;L,a;L',a') 



E 

ill 

x {[N LO G* Ma ^ LI {^-E + hw L o) 
+ (N LO + l)G* L . a , tL , (q,| ; E - fiw LO )] 

+ \[G<L: a 'M(<lW:E-f^L0) 

-G< L . a , !L ,(<i f E + fiw L o)}} (31) 

where N^q is the Bosc-Einstcin distribution for equilib- 
rium bosons with energy E p h on = Hujlo an d at lattice 
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temperature To. M pop is a basis dependent function of 
the coupling parameters, spatial structure and momen- 
tum transfer. 

For low energy (elastic) scattering with acoustic 
phonons and high lattice temperature, the expression 
for the equilibrium phonon propagator can be simplified 
to provide the (block-)diagonal, momentum independent 
self energies 

k 

A detailed derivation of the electron-phonon self energies 
for zinc-blende structures can be found e.g. in ((HI). 



C. Quantum kinetic equations 

Within the NEGF formalism, the steady state equa- 
tions of motion for the Green's functions are given (in 
matrix notation) by the Dyson's equations 

G R (k, 25) = [{G R (k, E))- 1 - £ R (k, E) - £* B (k, 25)] , 

(33) 

G*(k, E) = [{E + i v )l - Hoik)}- 1 , (34) 

G^(k, 25) = G R (k, E) (s^(k, E) + S^ B (k, 25)) G A (k, E), 

(35) 

G A (k,E) = (G R (k,E))l (36) 

Together with the expressions for the self energies from 
boundaries and interactions, and the macroscopic Pois- 
son equation 



layer L is 



eo 



dz 



e( Z )-U( Z ) 



= n(z)-p(z)-N dop (z), (37) 



relating the Hartree potential U(z) to doping density 
Ndop(z) and the carrier densities derived from the Green's 
functions, these form a closed set of equations for the 
latter that have to be solved self-consistently. To lower 
the computational costs, the recursive Green's function 
mctho d 20 i 25 ' 26 is applied, and only the first (block)off- 
diagonal of the self energies is considered. 



<^ = --XKY<j §tr{G<l\k;E)}. (40) 

k 

where A denotes the cross section area and A the layer 
thickness. The current density passing from layers L to 
L + 1 is 



J 



n(p) 



2e ^ f dE 
HA 



tL+^LGf^l^E)}. 



E. Absorption 



(41) 



The absorption of a given layer in a illuminated het- 
erostructure can be derived in terms of the microscopic 
interband polarization II C 



27 



a L (fuv) 



4tt 



-9{(Il£) £;L (g = (),&«;)}, (42) 



yfeAAchuj 
where 

9?{n*(k, 25)} = [H> (k, 25) - n< (k,25)] (43) 
= -- [H< (k, 25) -n<(k, E)\. (44) 

and 



n|(o, 25) = -2i J d ^L q \M2M\ 2 

xGf(d,E')G^(d,E' -E). (45) 

The incoming photon flux (/) 7 , after passing through lay- 
ers Li, 2/2) ••! Ln, is reduced by the absorptivity 

N 

a 7 = 1 - CXp[- E a £n(^7)A]; </>j,abs = 4>i a j (46) 



where </> 7 , a f)s is the absorbed photon flux. 



F. Computational scheme 



D. Carrier and current density 

The local density of states (LDOS) at layer L is given 

by 

p L {E)=Y,tr{A L , L {k ] E)}, (38) 

k 

A = i(G R — G A ), (39) 

where A is the spectral function and the trace is over 
orbital indices. The averaged electron (hole) density at 



After choosing an initial potential profile (e.g. from 
the depletion approximation) , the boundary self energies 
are calculated and used in the Dyson equation ([33| for 
the retarded Green's function G R , followed by the eval- 
uation of the Keldysh equation (f35|) for the correlation 
functions G^. These Green's functions provide an up- 
date of the scattering self energies T, R ^ (|28|) - (|32|) and the 
values of density and current. The new self energies are 
again used in the equations for the Green's functions, and 
this self-consistency iteration is continued until conver- 
gence is reached. Since the calculation of (photo)current 
is central to this work, its convergence is used as the 
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aborting condition instead of that of the Green's func- 
tions or self energies^. To obtain the built-in electric 
field, but also in cases where charging effects cannot be 
neglected (e.g. deep wells at large bias), Poisson's equa- 
tion is solved in an additional self-consistency loop using 
the densities from the NEGF and providing an update to 
the Hartree potential in the ballistic Hamiltonian. The 
computational scheme is represented in Fig. 



is 1000 W/m 2 (~ 1 sun) and the cross section is A = 1 
cm 2 . The photon energies are chosen in the range of the 
confinement level separation between the two band gap 
values, such that the contact and lead regions are non- 
absorbing. 



A. Local density of states 




FIG. 2: Computational scheme for the calculation of physical 
quantities from Green's functions and self energies. The inner 
self-consistency loop connects the equations for the Green's 
functions and the self energies, while the outer loop provides 
the update of the Hartree potential from the solution of Pois- 
son's equation. 



Since the system is open, there arc no true bound 
states, and the formalism considers only states contribut- 
ing to current, i.e. connected to extended states with fi- 
nite amplitude in the contacts. Fig. shows the local 
density of states (LDOS) for a 25 ML well at k = 0. In 
this case, two sharp confinement levels are present and 
contribute to the photocurrent. The high lying state is 
only weakly bound and broadened, corresponding to a 
faster carrier escape as compared to the more strongly 
bound and sharper low lying state. In the case of strong 
scattering, phonon satellite peaks form next to the con- 
finement level peaks, as visible in the cut of the LDOS 
through the center of the well (Fig. [3Jd). In addition 
to the confined states, there is a variety of quasi-bound 
states and transmission resonances above the well, which 
influence the photovoltaic properties of the structure and 
might explain the enhanced absorption of QWSC ob- 
served at photon energies above the higher bandgap. One 
can further observe a kind of "notch" states between well 
and the corresponding contacts, as are usually observed 
in the presence of barriers. If scattering in the leads is 
neglected, a stripe type interference pattern forms due to 
reflection of carriers injected below the band edge at the 
contacts, above which the LDOS acquires the expected 
uniform value of the quasi-continuum, which however is 
still affected by the presence of the well. 



Optical transitions, absorption and 
photocurrent response 



III. RESULTS AND DISCUSSION 

The following results for generic single quantum well 
(SQW) pin-diodes were obtained using the two band sp z - 
Hamiltonian with parabolic and isotropic tranverse dis- 
persion discussed in the appendix. Table [J shows the set 
of microscopic parameters used in the simulations. To 
lower the computational burden, short structures of 70- 
100 nm with reduced energy gaps of 0.5 eV (well) and 
0.9 eV (barrier) are investigated. The band offsets of 
barrier and well material are chosen to resemble those of 
the GaAs-AhjGai-^As System with x ~ 0.3, i.e. 0.25 eV 
for the conduction band offset and 0.15 eV for valence 
band discontinuity. The contacts are made of 50 mono- 
layers (ML) of high bandgap material with strong doping 
(Nd, a = 10 18 cm" 3 ). Between contact and active device, 
intrinsic buffer regions of 60 ML are inserted. The calcu- 
lations are performed at 300 K, the illumination intensity 



The different optical transitions between confined 
states, quasi-bound states, higher resonances and the 
continuum can be identified in the photocurrent response 
(Fig. [3h), which at short circuit conditions corresponds 
to the external quantum efficiency, i.e. the short circuit 
current normalized by the incoming photon flux. Since 



TABLE I: Material parameters used in simulations 





barrier 


well 






0.75 


0.5 


s-orbital onsite energy 


E p 


-0.15 





Pz-orbital onsite energy 




2.8 


2.5 


layer coupling element 


m* Cb /m 


0.1087 


0.067 


effective mass in conduction band 


my b /mo 


0.29 


0.23 


effective mass in valence band 


c 


12.2 


13.1 


dielectric constant 


M 


1 


1 


magnetic permeability 
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FIG. 3: (color online) a) Local density of states (LDOS) at k = for a 25 ML SQW pin-diode at Vwas = —0.01 V: quantum 
confinement leads to the formation of quasi-bound states and higher transmission resonances in the well region, in addition 
to the stripe- like interference pattern due to the built-in field (strong band bending); b) LDOS at the well center and optical 
transitions between confinement levels: the quasi-bound states near the well edge show the characteristic broadening associated 
with shorter carrier dwell time, as compared to the sharp deep and strongly-bound states; c) Photocurrent response (per) and 
absorptivity (abs): step like and square root like dependence on the photon energy below and above the higher band gap value, 
reflecting the density of the states participating in the corresponding optical transitions, i.e. confinement level to confinement 
level, confinement level to quasi-continuum and quasi-continuum transitions, respectively. 



the devices considered in this investigation are short, 
photocurrent is limited by the absorption, and it is there- 
fore essential to normalize physical quantities to the ab- 
sorptivity (Fig. [3b) in order to allow a comparison of 
different structures. 



C. Current spectrum and IV-characteristics 

There are two contributions to the total current in illu- 
minated QW pin-diodes: dark current, corresponding to 
the diode current driven by an applied external bias, and 
the photocurrent originating from the photogeneration 
of electron-hole pairs. Resolution in space and energy of 
the current in the QW region (Fig. [3t) allows the dis- 
tinction between the two components. The diode current 
occupies a narrow region above the band edge at the con- 
tacts, it is constant over the hole device and its spectrum 
reflects the density of states and the distribution of the 
carriers in the contact reservoirs from which they are in- 
jected, broadened by scattering with acoustic phonons, 
and relaxed towards lower energies by interaction with 
polar optical phonons. In the absence of interband re- 
combination, the diode current is conserved for electrons 
and holes separately. Photocurrent, on the other hand, 
is driven by the excitation of carriers from the opposite 
band, and current conservation^ thus holds only for the 
sum of electron and hole contributions, but not for the 
separate components, which increase towards the respec- 



tive contacts (Fig. \%jp) and differ also in their spectrum 
(Fig. (Us). The photocurrent spectrum reflects the joint 
density of states of the dominant transition between con- 
finement levels. Unlike the LDOS in the well, the current 
spectrum shows a strong asymmetry between electrons 
and holes: in the conduction band well, the main contri- 
bution to current comes from the higher level, while it is 
the lower one that dominates the current in the valence 
band well. This demonstrates the impact of carrier es- 
cape probability on the current, the latter no longer being 
characterized by the LDOS alone as in bulk structures. 

Fig. [5^, shows the current-voltage characteristics for 
the 25 ML SQW structure. Near short circuit conditions 
(Fig. Eb), current is purely photocurrent. At increasing 
bias, the diode current evolves exponentially (Fig. [5b,d), 
showing the specific spectrum of the injected carriers and 
the effects of scattering in terms of phonon satellite peaks 
towards the band edge. The spectrum of the photocur- 
rent is modified due to the Stark effect. 



IV. SUMMARY AND CONCLUSIONS 

We presented a microscopic model for the consis- 
tent description of generation and transport processes in 
semiconductor quantum well structures under monochro- 
matic illumination and in the radiative limit. Based on 
the NEGF formalism for a tight-binding Hamiltonian, it 
provides access to non-equilibrium phenomena in quan- 
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FIG. 4: (color online) a) Spatially resolved photocurrent spectrum (at zero bias voltage) in the QW region and b) at the 
interface to the n-contact (electrons) and to the p-contact (holes): the spectrum reflects the joint density of states for the 
contributing transitions between the confinement levels, modified by the probability for escape, which is suppressed in the 
case of the deep electronic level, c) Electron and hole components of the photocurrent grow towards the respective contacts, 
while the total current is conserved. In the dark, the bands are uncoupled, and current is conserved for the two carrier species 
separately. 
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FIG. 5: (color online) a) IV-characteristics for a 25 ML SQW 
structure and the current spectrum at the lead-device inter- 
face for b) V (short circuit conditions), c) -0.26 V (near 
the maximum power point), d) -0.32 V (near the open circuit 
voltage). The spectrum of the exponentially increasing diode 
current reflects the density of states and the distribution of 
the carriers in the bulk contacts, modified by the effects of re- 
laxation due to inelastic scattering in the active region, which 
leads to the formation of phonon satellites (weakly recogniz- 
able near the band edge). 



The insights into the photovoltaic performance of specifi- 
cally coupled multi-quantum- well structures, gained from 
the application of the presented approach, are the subject 
of current investigations and will be published elsewhere. 
Future work will also include a microscopic treatment 
of the main nonradiative recombination processes, which 
are Auger and trap recombination. For comparison with 
experiment, a more realistic band structure model will 
be used. In order to account for optical processes in ex- 
tended structures, such as photon recycling, the spatial 
variation of the light intensity needs to be considered, 
which can be accomplished by the solution of an addi- 
tional Dyson equation for the photon propagator con- 
taining the microscopic polarization function. Investiga- 
tions of hot carrier effects will require a corresponding 
treatment of the phonons in the quantum region. 
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APPENDIX A: TWO BAND TIGHT-BINDING 
MODEL 



turn confined structures subject to interactions and there- 
fore supports the investigation of the microscopic pro- 
cesses governing the physics of quantum well solar cells. 



The simplest tight-binding model to describe the con- 
duction and valence band structure of III-V semiconduc- 
tors like e.g. GaAs is the diatomic model with a two- 
orbital basis2£. In this model, the two-band dispersion 



8 



is reproduced approximately by placing an s-type or- 
bital on the cation (Ga) and a p z -type orbital on the 
anion (As). Fig. [5] shows a projection of the zinc-blende 
lattice onto the (001) direction, with the corresponding 
intra- and interlayer couplings U ac = U ca = V sp and 

V ac = V ct 

and E„ = 



= — V sp , and the "on-layer" energies E c = E s 
E p . For further simplification, the transverse 



A = a L /4 A = a L /2 



L - 1 



Uca 

L 



L + l 



FIG. 6: Tight binding elements for a zinc-blende lattice in 
(001) direction. 



band structure is approximated by an isotropic parabolic 
dispersion, with the electronic effective mass replaced by 
the hole effective mass in the case of the valence band. 
What remains in the direction of propagation is a sp z 
model with k = 0, which is equivalent to a linear chain 
of atoms, with the anion layers at z = na^ / '2 the cation 
layers at z = (na^/2 + ol/4), where n G N labels the 
monolayer of thickness A = cll/2, with the lattice 
constant of the binary compound (i.e. 5.65 A in the case 
of GaAs). In terms of the parameters introduced above, 
the elements of the TB-Hamiltonian 



/ 



H TB — 



\ 



tL-lL-2 DL-ilb tL-lL 

thL-1 D^tb tLL+l 

tL + lL Dz+llb ^L+lL+2 



V 



including parabolic transverse energy, are given by 



/ 

(Al) 



Dl = H nn — 



2m*. 



-v 



VsT) E v 



sp 

h 2 k 2 



(A2) 



N 



tLL+l = H nn+ l - ^ q j , tLL-1 = - ^ ^ P 

(A3) 

where rn* el i hl is the effective electron and hole mass, re- 
spectively. This yields the bulk Hamiltonian 



H(k, k z ) 



E, 



h A k A 



2iV sp sm{k z ^) 



-2iV sp sm(k z ^) 



En 



(A4) 



which for k = gives rise to the dispersion relation 



det[# (fc,) - E] = ^E(k z ) = ^ 



E p + E s 



±J(Ep-E s f + 16V s %sm 2 (k z ^ 



(A5) 



For the integration over transverse momentum, the 
isotropic one dimensional approximation 



E 



A 



BZ„ 



d 2 k pa ^- j dkk, (A6) 



is used, where A is the device cross section, BZ\\ is the 
projected Brillouin zone and k = |k| is the absolute value 
of the transverse momentum. 

The tight-binding parameters are related to the longi- 
tudinal effective mass m* through the longitudinal dis- 

l 

, with E z {k z ) re- 



persion relation, as m* 
suiting from the secular equation 



Tr 
m 



ilk' 2 



det[H ± (k z ) - E z )] - 0, H ± %) = (k, k z 



21,2 



K 2 k 

2^* 



(A7) 

From Eq. (|A5p . one finds the relation between the effec- 
tive mass at the T-point and the coupling element V sp , 



a 2 V s 2 p q 2 



n 2 



m \2\E S -E p \q 



V sp 



2E n 



where E g = \E S — E p \ is the energy gap. 



a V to* 1 tooq 

(A8) 



APPENDIX B: BOUNDARY SELF ENERGIES 
FOR MULTIBAND TIGHT-BINDING 
TRANSPORT MODELS 

To properly model the effect of semi-infinite bulk at 
the lead-device interface, the interface Green's function 
has to be linked to the propagating and evanescent states 
in the leads. The total electron wave function expressed 
in terms of the Bloch sum of the anion (a) and cation (c) 
states as a linear combination of planar orbitals \a, L, k) 
is given by Eq. |3|). In the planar orbital basis, project- 
ing onto the atomic orbitals a' located at layer L, the 
Schrodingcr equation for the contact Bloch states reads 



^2(a',L,k\H\a,k z )=0, 



(a',L,k\H\a,k z ) = (a',L,k\H\a,k z ) 



— E(a , L, k\a, k z 



(Bl) 



(B2) 



For a tight-binding Hamiltonian coupling m neighboring 
layers, which is of the form 



H(k,k z )= J2 H°{k)e l ° k : 



A 



(B3) 
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where H rT (k) represents a matrix which couples a given 
layer to the cr-th neighboring layer and A is the layer 
spacing, and defining 



C^s^X a = -m,..,m, 
Eq. ()B2[) can be written as 



(B4) 



J2 H'cr + ir 



0, (B5) 



where it was used that C™ = e ifc * A C m_1 . 
For a nearest neighbor Hamiltonian (m = 1). the pro- 
jected Schrodinger equation is recast into 



a—lr-ta—1 



H"~ C 
which, using C 



iJ^C" 7 + iJ CT+1 C" T+1 = 



(B6) 



crrtl 



_ e ±ifc.A^<o-^ can wr i-tten as 
/PC" 7 + i^+V^C" 7 = 0. (B7) 



This equation can then be transformed into an eigenequa- 
tion for the propagation factors A = e lk * A and the lead 
Bloch states in local orbital basis: 



TC L = XCr 



C 



L+l 



(B8) 



with Cl = ( ^ a j and T = T c T a , where T a and T c are 
the atomic layer transfer matrices defined as 



Th = 



(b) 

1,1-1 



H 



H, 



(b) 
1,1-1 



H. 



(b) 
1,1+1 



(b = a, c) 



(B9) 



with the matrix elements given by (I denotes the atomic 
layer) 



(b) 

l,l—l,a.Oi' 
(b) 

l.l.a.a' 



h. 



(b) 

l,l+l,a,a' 



(a,l,k\H\a',l-l,k), 
(a,l,k\H\a',l,k) - E6 
(a,l,k\H\a',l + l,k). 



(BIO) 
(BH) 
(B12) 



The eigenstatcs \ and eigenvalues A 



ik,A 



ofEq. (TB8 



correspond to the bulk modes propagating (real k z ) or 
decaying (complex k z ) to the left (3?(fc z ) < 0) and to the 
right ($t(k z ) > 0), respectively. For an iV^-band model 
with a two atom basis, there are states Xv propa- 

gating or decaying to the right (y = +) and to the left 
{v = — ), respectively. At a given layer L, the compo- 
nents for left- and right travelling waves can be written 
as 



Cl± = U±C±, 



(B13) 



where C„ is a vector containing the expansion coeffi- 
cients, and 







Ac)N b /2 



(B14) 



The corresponding expression for the adjacent layer L + l 
is 



Cl+i± — U±\f 1 C±, 
with the propagation matrix 



/ e lk ' A 



(B15) 



A, 



oifcf 6/2 A 







V 



ikiA 







(B16) 



e^ b/2A / 

The relation between the two layers follows as 

C(l+i)± = F±C L± (B17) 

with 

F± = U±Af 1 U ± 1 . (B18) 

Relation (|B17|1 can be used to derive the retarded Green's 
function g rR at the right boundary (L = 1) of the un- 
coupled semi infinite left lead, i.e. for the case where all 
the couplings to the right are set to zero. The equation 



[(E + i V )l-H ]g rR = l 
yields for the boundary element 

[(E + i V )l-Di}gl*-ti.,ig r o *! = 0. 
Eq. (|B17j) provides the relation 

rR _ W -lrR 



(B19) 
(B20) 
(B21) 



which determines the left lead boundary Green's function 
in terms of the bulk modes as 

g r 1 R = [Et-D 1 -ti ;0 FI 1 ]- 1 (B22) 
= [El -D 1 -E^f] -1 , (B23) 

and providing thus an expression for the (left) retarded 
boundary self energy S^f . 
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